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ABSTRACT 

We present a new analysis of the light curve for the secondary star in the eclipsing binary 
millisecond pulsar system PSR B 1957+20. Combining previous data and new data points at 
minimum from the Hubble Space Telescope, we have 100% coverage in the R-band. We also 
have a number of new Kg -band data points, which we use to constrain the infrared magnitude 
of the system. We model this with the Eclipsing Light Curve code (ELC). From the modelling 
with the ELC code we obtain colour information about the secondary at minimum light in 

o 

BVRI and K. For our best fit model we are able to constrain the system inclination to 65 ± 

o 

2 for pulsar masses ranging from L3 - 1.9 Mq. The pulsar mass is unconstrained. We also 
find that the secondary star is not filling its Roche lobe. The temperature of the un-irradiated 
side of the companion is in agreement with previous estimates and we find that the observed 
temperature gradient across the secondary star is physically sustainable. 
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1 INTRODUCTION 



The 



bin ary 



millisecond 



pulsar 



(MSP) 
the 



PSR 



However, the means b y which such a scenario could occ ur is still 
the subject of debate lEichler & Levinson 1988 & 1991; see also 



B 1957+20 jprucht er. Stinebri ng & Taylor lysai) is tne origi- 
nal and one of the best studied members of its class. It consists of 
a 1.6 ms radio pulsar orbiting a companion of mass no less than 
O.O22M0 , in a binary of orbital period 9.17 hours. For 10% of this 
orbit, the radio emission from the pulsar is eclipsed: the eclipsing 
region is considerably larger than the Roche lobe of the companion 
star, suggesting a wind of material from the secondary star, due t o 
ablation by the impinging pulsar radiation jpruchter et al. 1988h . 
The optical counterpart to PSR B 195 7+20 was discovered by 
iKulk arni. Diorgovski & Fruchter (1988): subsequent observations 
by ICallanan, van Paradijs & Rengelink ( 1995) found the optical 
counterpart to vary by a factor of 30-40 in flux over the course of 
the orbital period. These observations are hampered by the close 
proximity of a line of sight "contaminator", only ^ 0.7" away. 

At the time of its discovery it was assumed that PSR 
B 1957+20 was the missing link between low-mass X-ray binaries 
(LMXBs) and isolated millisecond pulsars. It was suggested that 
the high energy pulsar radiation could evaporate the companion 
star (Ruderman et al. 1989a & 1989b) leavin g behind an isolate d 
millisecond pulsar like PSR B 1937+21 feacker et al. 1982h . 



iBhattacharva & van den Heuvel 199lL iPhinnev & Kulkarni 19941 
and i Lorimer 20011 for thorough reviews of the formation and 
evolution of MSPs, binary MSPs and the current status of MSP 
research respectively). 

PSR B 1957+20 is a member of a class of binary pulsar sys - 
tems, the Black Widow Pulsars (e.g. iKing. Davies & Beer 20031) . 
These are systems with secondaries of mass typically less than 
O.O5M0 and orbital period less than 10 hours. To date 15 such 
systems have been identifiecQ with 13 of these residing in globular 
clusters and only 2 situated in the field. Radio eclipses have been 
detected in approximately half of the cluster systems and in both of 
the field systems. Optical light curves exist for only the field sys- 
tems, PSR B 1957+20 (Callanan ~al. 19951) and PSR J205 1-0827 
( iStappers et al 19991 20011) . 

Previous attempts at modelling the optical light curve of PSR 
B 1957+20 had one major limitation: the companion was unde- 
tectable at minimum. As such it was impossible to tightly constrain 
important system parameters such as the inclination or the degree 
of Roche lobe filling. However, a number of R & I-band images 
of the optical counterpart at minimum were subsequently obtained 
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by the Hubble Space Telescope (HST). Furthermore, we have re- 
cently acquired a number of Ks-band images of the near infrared 
counterpart. In this paper we combine these data for the first time to 
model the light curve to obtain tight constraints on the inclination 
and Roche lobe filling fraction. 



2 DATA 

2.1 Optical Photometry 



Table 1. Ks magnitudes of PSR B 1957+20. 



4' 


Secondary Ka mag 


Contaminator Ks mag 


0.5 


17.8 ±0.1 


18.06 ±0.06 


0.764 


18.20 ±0.05 


17.98 ± 0.04 


0.774 


18.27 ±0.10 


18.13 ±0.08 


0.788 


18.35 ±0.06 


18.02 ± 0.05 


0.824 


18.73 ±0.09 


18.11 ±0.07 



The data set consists of B, V and R-band images taken with the 
William Herschel Telescope ( WHT) at La Palma on the nights of 
1989 July 2-5 (see Callanan et al. 19951 and references therein for 
details of the observations), along with two pairs of R (1994 Aug 
30, Oct 28) & I-band (1994 Sep 03, Oct 27) data points taken at 
minimum with the HST. 

We re-reduced the WHT data using the IRApH implementa- 
tion of DAOPHOT (.Stetson 1987h. We und ertook this task as, in their 
original analysis, ICallanan et al. (19951) used only approximately 
half of their data due to poor seeing during 2 of the 4 nights of their 
observing run. Our re-analysis allows us to use an additional night 
of data. In the end we had 41 useful R-band data-points covering ~ 
85% of the orbit, 39 V-band data points covering ~ 70% of orbital 
phase and 7 useful B-band points taken near maximum. 

Two R-band & 2 I-band data points were also obtained during 
eclipse with the Wide Field Planetary Camera 2 using the F675W 
and F814W filters respectively. The exposure time was 600s in each 
case. These images were corrected for cosmic ray hits and the ob- 
ject magnitudes were calculated using the QPHOT task in IRAF. 
The F675W and F814W magnitudes were converted to R and I- 
band magnitudesjii the J ohns on system following the prescription 
given bv lHoltzman et al. (19951) . 

The WHT and HST data was then phased ac cord- 
ing to the radio ephem erides of lRvba&Tavlor(I99lh and 
lArzoumanian et al. (1 9941) respectively. When combined, these 
points provide the first complete optical light curve of PSR 
B 1957+20, in addition to colour information throughout the orbital 
cycle. 
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Figure 1. Ks band image of the PSR B 1957+20 system. The arrow indi- 
cates the position of the secondary with the line of sight contaminator (0.7" 
sepai'ation) lying to the north east. The exposure time was 180s. 



2.2 IR photometry 

Our observations consist of a series of Ks-band images obtained 
with the 6.5m Magellan Baade telescope at Las Campanas Obser- 
vatory on 2004 July 23 and 2005 September 13 using the PANIC 
camera. These were dark-current subtracted, flat fielded, back- 
ground subtracted and combined using standard IRAF routines. 
In total we obtained 4 images totalling approximately 18 min- 
utes on source exposure time. We display one of our K^-band 
(1.99/im ^ A ^ 2.30/im ) images in Fig. [T] The companion 
star to the pulsar is easily resolved. 

Photometry was carried out in the same manner as the WHT 
op tical data. The image was calibrated using the standard star P565- 
C jPersson et al. 1998h and cross checked by comparing a num- 
ber of stars in our field with those in the 2MASS catalogue. Th e 
data was phased using the ephemeris of lArzoumanian et al. (1994h . 

^ IRAF is distributed by the National Optical Astronomy Observatories, 
which are operated by the Association of Universities for Research in As- 
tronomy Inc., under cooperative agreement with the National Science Foun- 
dation. 



which at the time of our observations was accurate to at least 1 
second. Our final K^-band photometry is displayed in Table[T] 

Previous attempts at IR photometry of this system were un- 
able to resolve the pulsar from the line of sight contaminator and 
as such only the K -band magnitude of the unresolved combination 
was obtained (see lEales et al. 19 9C/). However, our new observa- 
tions allow us to subtract t he magnitude of th e contaminator from 
the combined mag nitude of lEales et al. (199(]|) . yielding a K^-band 
magnitude of the system at maximum of 17.8 ±0.1. 



3 THE ELC MODEL 

To m odel these data w e used the ELC light curve modelling 
code dOrosz et al. 2OO0h . The ELC code is ideally suited for this 
type of system as it incorporates the NEXTGEN low tempera- 
ture model atmosphere tables, which are critical for systems like 
PS R B 1957+20, with a companion of likely temperature ~ 3000 
K jpruchter et al. 19951) . The ELC code also allows one to fit light 
curves on a one by one or simultaneous basis. In our case this al- 
lowed the fitting of the BVRI and K-band light curves simultane- 
ously. 
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Table 2. Orbital Parameters of PSR B 1957+20 system used in the ELC 
modelling. 



Parameter 


Value 


Companion Mass^ 


0.022 Mq 


Orbital Period^ 


9.17 hr 


1^§ ^spindown 


35.20 erg s-l 


Companion Effective Temp'^ 


2800 K 


Inclination'^ 


50 - 80° 


f 


~ full 



1 Fruchter et al.l988 ; ^ Fruchter et al. 1995 
^ Callanan et al. 1995 




3.1 The Model 

The ELC program requires a number of input parameters be- 
fore modelling the light curve: the initial parameters used are 
given in Table |2] The temperat ure dependent gravity darken- 
ing exponents of Iciaret (200d) were used. We initially at- 
tempted to model the system as a blackbody (T ~ 2800 
K); however, the results were unsatisfactory. While the code 
had no problem in fitting the observed light curve at maxi- 
mum, it was completely incapable of reproducing the observed 
minimum (the model was consistently too luminous during 
eclipse). We then employed the NE XTG EN model atmospheres o f 
iHauschildt, Allard & Baron (1999ah & iHauschildt et al. (1999bh . 
and using our blackbody model as our starting point, we proceeded 
to model the light curve. The pulsar mass was set to the canon- 
ical value of 1.4 Mq. We then varied the following parameters: 
inclination and mass ratio of the system, Roche lobe filling fraction 
(/), temperature and bolometric albedo (a) of the secondary star 
and the irradiating luminosity. The geneticELC algorithm (based 
on the PIKAIA routine of Charbonn eau 1995) was used to search 
for the best fit values. The best fit R-band model is displayed in 
Figure |2] We see that there is excellent agreement between the fit 
and the data (xt = 1.06). The largest deviations occur at orbital 
phases > 0.65, but even these are well within the errors. This dis- 
crepancy is due to the relative faintness of the companion at these 
phases, and poorer seeing conditions during these observations. As 
a check on the validity of the model we used our limited K^-band 
data, as displayed in Figure |3] We see that the fit agrees with these 
data very well. 

Given that the mass of the pulsar is currently unknown, al- 
though most evolutionary scenarios suggest that it will have ac- 
creted a few 1/10*'' of a solar mass from the secondary, we decided 
to repeat the above procedure for a number of other primary masses 
in the range 1.3 < Mmsp < 1.9 Mq, to investigate the effect of 
the pulsar's mass on our estimates of the mass ratio and orbital in- 
clination of the system. 



4 RESULTS 
4.1 Inclination 

At a given pulsar mass the inclination was constrained to within 



± 1.2 , i.e. for a pulsar of mass 1.4 Mq, ; = 64.4 



+ 1.3 



(see 



Figure|4l(, and overall for the above range of pulsar masses we find 

o o 

the inclination of the system to be in the range, 63 ^ z ^ 67 , at 
the 3a level. 



Figure 2. The best fit to the combined R-band data with residuals (inset). 
Two orbital phases are displayed for added clarity. The pulsar mass is 1 .4 
Mq . The best fit incfination is / = 



■-64A+\ l (3(7) with a x5 = 1.06 



18.5 
19 




Figure 3. The simultaneous fit to the K-band light curve corresponding to 
the R-band fit in Figure|2] 



4.2 Pulsar Mass 

The value of xi exhibited only a nominal increase as the mass of 
the pulsar was increased from 1.3 - 1.9 Mq: hence our models are 
unable to constrain this parameter. 



4.3 Roche lobe filling factor 

At no point in our attempts to model this system were we able to 
obtain an acceptable fit for a secondary filling its Roche lobe. For 
our models using the NEXTGEN model atmospheres the value of 
/ was approximately constant, 0.81 ^ / ^ 0.87 (3a level), as we 
varied the mass of the pulsar between 1.3 and 1.9 Mq. Hence the 
secondary is tightly constrained as not currently filling its Roche 
lobe. 



4.4 Temperature of the Secondary at maximum and 
minimum 

We obtained a value ofT = 2900 ± 1 10 K (3a level), for the effec- 
tive temperature of the un-illuminated side of the companion star 
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Figure 4. The graph of vs (' for a pulsar of mass 1.4 Mq. Minimum 
occurs for xt = 1-06 and an inclination, i ~ 64.4 °. The 68%, 95% & 99% 
confidence levels are illustrated. 

for pulsar masses in the range 1.3 - 1.9 Mq. For individual pulsar 
masses the 3(t error was ± 90 K i.e. for a pulsar of mass 1.4 Mq 
an effective temperature of T = 2900 ± 90 K was obtained. The 
corresponding temperature at maximum is T = 8300 ± 200 K (3cr). 

From our modelling, we have obtained the magnitude of the 
secondary at maximum (in I) and during eclipse (in B, V and K); 
see Table [3] These provide us with colour information about both 
the cool side and the irradiated side of the pulsar's companion. 

We find that at least ~ 70% of the spin-down energy of the 
pulsar is required to produce the observed heating effect and that 
this percentage is independent of the mass of the pulsar. The bolo- 
metric albedo of the system was found to remain close to a value 
of 0.5 for all our models, which ensured that the secondary was 
convective. 



5 DISCUSSION 

We have modelled the light curve of the PSR B 1957+20 system 
with the ELC code and we find the system to be accurately mod- 
elled by a highly irradiated secondary. The inclination of the system 

o o 

is measured to be 65 ±2 for a pulsar in the mass range 1.3 - 1.9 

Mq. 

The optical/IR lightcurves display no evidence for the pres- 
ence of a contribution from the intra-binary shock. If such a shock 
did contribute in a non-negligible manner to the optical/IR flux 
from the system, we would expect to observe this in the form of an 
asymmetrical distortion of the lightcurves, which is not observed. 
In contrast the highly symmetrical nature of the lightcurves is strik- 
ing evidence that the modulation is the result of the emission from 
the heated face of the secondary star. Bogdanov et al. (2005) ob- 
served variable non-thermal X-ray emission attributed to the intra- 
binary shock in the binary millisecond pulsar 47 Tuc however, 
an extrapolation of this emission to optical wavelengths demon- 
strated that it contributed negligibly here. Recent XMM-Newton 
observations (Huang et al. 2007) have tentatively detected similar 
emission in the PSR B 1957+20 system, although our observations 

3 also known as PSR J0024-7204 



show it to have an insignificant contribution in the optical, as ob- 
served in 47 Tuc W. 

5.1 The Roche lobe filling factor 

The Roche lobe filling factor is constrained to be greater than 
80%, for our best fit model. This result is in agreement with 
past estimates that r equired the seconda r y to be close to fill- 
ing its Roche lobe ^Aldcroft et al. 1992*; Applegat e et al. 19941 : 
iBrookshaw et al. 1995 : Callanan et al. 1995). Previo us es timates of 
the mass loss rate in this system jFruchter & Goss 19921) required 
the system to be close to filling its Roche lobe, as the measured 
density of the eclipsing material was too tenuous to account for 
significant mass loss. They claimed that if the secondary neared its 
Roche lobe, material could easily leave the stellar surface and re- 
main in the orbital plane; this would explain the low density of the 
observed material. 

5.2 Temperature and albedo of the secondary star 

The effective temperature of the un-illuminated side of the sec- 
ondary, T = 2900 ± 1 10 K, is in excell ent agreement w ith the pre- 
vious estimate of T = 2800 ± 150 K jFruchter et ai~1995.) and is 
corroborated by its agreement with the R - 1 colour temperature ob- 
tained via the HST (~ 3000 K. ICox 2000l) . The temperature derived 
from the colour information at maximum is Tmax = 80001Jqqq K. 
This compares with a temperature of 8300 ± 200 K from the ELC 
models. The error in the colour temperature at maximum is domi- 
nated by the large uncertainty in the extinction in the direction of 
PSR B 1957+20. 

The bolometric albedo of the secondary, d efined as the r atio of 
the reradiated energy to the irradiance energy (^ Wilson 1990 ^, was 
found to favour the convective case (a ~ .5 ) as one woul d ex- 
pect for a secondary of such a small size ( iKing et al. 20o3) . We 
note here the modelling of the analogous system PSR J205 1-0827 
( Stappers et al. 2001), in which the percentage of the pulsar's spin- 
down luminosity, which is re-radiated in the optical by the sec- 
ondary, was determined to be in the region of 30% - 45%. This is 
consistent with the the incident spin-down luminosity that our mod- 
els require for PSR B 1957+20 (^ 70%), given the above albedo. 

As a check, a number of models were constructed in which the 
bolometric albedo was set to 1.0 (radiative secondary). In this case, 
the best fit value of the irradiating luminosity is found to be lower 
as expected but the associated temperature of the cool side of the 
secondary is inconsistent with the observed colours. 

5.3 Temperature Gradient 

It is clear from the light curve of PSR B 1957+20 that a large tem- 
perature gradient is required between the heated and cool hemi- 
spheres of the companion star. To test if this is physically sustain- 
able, we decided to model the heat flow along the surface of the 
secondary in more detail. 

A two-dimensional model of the irradiation of PSR B 1957+20 
was simulated using a modifi ed version of the code described in 
iBeer & Podsiadlowski (200 2j . The code uses a polytropic equa- 
tion of state and only hydro-dynamical effects were initially con- 
sidered. The irradiation causes a stress on the stellar surface which 
drives a sub-sonic circulation. Once this circulation pattern was 
found thermodynamic effects were considered. Matter in the di- 
rectly illuminated region was heated and the advection of this mat- 
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Table 3. Optical and IR magnitudes of the companion to PSR B 1957+20 at maximum and during eclipse. 





B 


V 


R 


1 


K 


Max 
Min 


21.08 ±0.05 
28.1 ± 0.1 


20.16 ±0.05 
26.2 ± 0.1 


19.53 ± 0.05 
24.6 ± 0.1 


18.79 ± O.OS 

22.52 ± 0.05 


17.8 ±0.1 
19.S ± 0.1 



The values in bold are predicted by the ELC models, all other values are measured directly from the photometry. 



ter was followed across the stellar surface. As the matter flowed it 
was allowed to cool radiatively. The resulting temperature distribu- 
tion was evolved in time until a steady-state solution was achieved. 
It was found that the heated material extended beyond the directly 
irradiated region but that not all of the un-illuminated portion of 
the star was heated. Consequently a large temperature gradient be- 
tween the illuminated and un-illuminated sides was found to exist. 

It may not be obvious how such a large temperature gradient 
can exist across the surface of the secondary. In fact the converse 
has also been argued. If the radiative cooling timescale is short 
then little or no energy would be redistributed jPahab 1974) . The 
reason why a large temperature gradient can exist is because the 
star is perturbed from hydrostatic equilibrium by the irradiation in - 
duced circulation as first noted by Kippenhah n & Thomas (igTsh . 
The circulation pattern attempts to distribute the energy due to heat- 
ing across the surface of the secondary. However, the circulation 
itself, produces inertia terms in the equations of motion which per- 
turb the star from hydrostatic equilibrium i.e. the pressure gradient 
is no longer in the same direction as the potential gradient. This in 
turn is what allows the large temperature gradient to exist across 
the surface, even of a star as small as the secondary present in this 
system. This effect is independent of whether the object in question 
is degenerate or non-degenerate. 



of its Roche lobe. In contrast a 50 Myr old 0.025 Mq brown dwarf 
would be expected to have a temperature of 2200 K and a radiu s 
approximately half the size of that observed dChabrier et al. 2000l) . 
[Applegate et al. (1994) have previously proposed a model in which 
the secondary star is heated to this temperature through tidal heat- 
ing; this model also has the advantage of naturally explaining th e 
orbital period variability observed by lArzoumanian et al. (19941) . 
However, this model requires the secondary star to be close to fill- 
ing its Roche lobe and given that we observe the secondary to be 
underfilling its Roche lobe by up to 20%, it is questionable if tidal 
heating would be an efficient heating mechanism in this scenario. 

The accreting millisecond pulsar SAX J1808. 4-3658 is also 
observe d to have a bloated, low mass ( '-^ 0.05Mq) companion. In 
this case lSildsten & Chakrabartv (2001 ) suggest that the secondary 
star is 'pumped up' to the bloated higher entropy state by the heat- 
ing effect of the thermal radiation emitted by the neutron star in qui- 
escence. One could envisage the secondary star in PSR B 1957-1-20 
being affected in a similar manner but in this case the heating would 
be caused by the incident spin-down radiation from the radio pul- 
sar. In reality the situation is likely to be a complicated interplay 
between both mechanisms which combine to produce the abnor- 
mal secondary present in this system. 



5.4 Nature of the secondary 

We can use the above colour/temperature information to constrain 
the nature of the secondary star. The colour information appears to 
rule out the possibility that the secondary is a white dwa rf. In their 
study of ultracool white dwarfs (T < 4000 K) Gates et al. (20041) 
found that white dwarfs at this low temperature typically have 
R-I < 0.5, whereas we find (R-I)o = 1.8 ± 0.3. The redden- 
ing in the direction of PSR B 1957-1-2 was calculated usin g the 
hydrogen column density estimate of 'Stappe rs et al. (2003h . Nh 
= (1.8 ± 0.7) X 10^^ cm in combination with the extinction 
curve o f ISavageetal. (19791) . Furthermore, if the secondary in PSR 
B1957-I-20 was a ~ 0.025 Mq white dwarf, one would expect a ra- 
dius of ~ 0.1 Rq - again in contrast with the value of ~ 0.3 Rq 
which we have de termined from our modelling. 

According to lBessell (199 li) R - 1 is the most reliable spectral 
type indicator for late M-type stars and using their diagrams of both 
spectral type & temperature vs R - 1 for a sample of late M dwarfs, 
we find that temperatures of between 2900 - 3100 K and spectral 
types of M4 - M7 are in agreement with our observed R-I. Hence 
the secondary appears to exhibit the colours of a late M-type dwarf, 
although a main sequence companion (Mcomp > 0.08 Mq) is ruled 
out on the basis of the mass function combined with our inclination 
estimate above. It is clear that the current mass of the secondary 
is well below the hydrogen burning limit of 0.08 Mq. Hence, the 
most likely current state of the secondary is that of a brown dwarf. 

The low mass secondary in this system has been observed to 
have a temperature of ~ 2900 K and a radius encompassing ~ 80% 



5.5 Comparison with other MSPs 

At this point we should also compare ou r results with the othe r 
ablating field system PSR J205 1-0827 ( Sta ppers et al. 199911 200 ll) . 
which is remarkably similar to PSR B1957-H20. The secondary in 
this system appears to be similar to that in PSR B 1957-F20. It has a 
mass of ~ 0.025 Mq and the temperature of the cool side has been 
measured to be ~ 3000 K. This system has also been found to be 
under-filling its Roche lobe, in this case by ~ 50%. Even though the 
orbital period of this system is only 2.4 hrs (in comparison to 9. 1 hrs 
for PSR B 1957-1-20), it is clear that a similar ablation mechanism is 
at work. 

Interest in these ablating systems has increased with the 
discovery of the accretion powered X-ray millisecond pulsars 
(AXMPs). There is evidence that the pulsar in these sys- 
tems is in the process of 'turn ing-on' as a radio pulsar. In 
the AXMP SAX J1808.4-3658 din't Zand et al. 1998h . consist- 
ing of a neutron star and a secondary of mass ~ 0.05 Mq 
teildsten & Chakrabartv 200 ll) . a similar process could be taking 
place. Once again the secondary is being heated, but in this case 
the c ool side of the seconda ry is only 1000 K cooler than the warm 
side. ICampana et al. (2004h have interpreted this as evidence that 
the companion star is being irradiated by the spin down luminos- 
ity of the pulsar. Further observat ions of quiescent AX MPs are re- 
quired to test this hypothesis (e.g. ICallanan et al. 20071) . 
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6 CONCLUSIONS 

The main aim of tliis paper was to constrain the inclination 
of the system as a precursor to a campaign of phase resolved 
spectroscopy, with the aim of measuring the mass of the pulsar. 
We have determined the inclination to within ±2° for a pulsar 
mass in the range 1.3 - 1.9 Mq. This should ensure that any mass 
determination will be limited only by the accuracy of the radial 
velocity measurements. We have shown that the temperature of 
the secondary agrees with previous estimates and the observed 
temperature gradient is physically sustainable. We also find the 
secondary to be under-filling its Roche lobe by up to 20%. 
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